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Abstract 

The usual e, <5-definition of the limit of a function (whether pre- 
sented at a rigorous or an intuitive level) requires a "candidate L" 
for the limit value. Thus, we have to start our first calculus course 
with "guessing" instead of "calculating" . In this paper we criticize the 
method of using calculators for the purpose of selecting candidates for 
L. We suggest an alternative: a working formula for calculating the 
limit value L of a real function in terms of infinitesimals. Our for- 
mula, if considered as a definition of limit, is equivalent to the usual 
e, 5-definition but does not involve a candidate L for the limit value. 
As a result, the Calculus becomes to "calculate" again as it was orig- 
inally designed to do. 

1 Introduction 

Let / : X — > K be a real function, where X C R and r G R be a non-trivial 
adherent (accumulation) point of X. Recall that a real number L is called 
the limit of / as x tends to r, in symbols, lim x ^ r f(x) = L, if 

(Ve G R+)(35 G M+)Vx G X)(0 < \x - r\ < 5 =^ \f(x) - L\ < e). 

This is the so-called e, 5-definition of limit. It is sometimes attributed to 
Cauchy but it appears historically for the first time in John Wallis's "Arith- 
metica Infmitorum" (The Arithmetic of Infinites) in 1655. In the next 250 
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years following Wallis's work, this definition was rejected and rediscovered 
many times until it was finally accepted by the mathematical community 
in the beginning of the 20th century. We should mention that at that time 
most of the results in Calculus were already discovered through infinitesi- 
mals. Contemporary mathematicians might be puzzled by the fact that it 
took so long for the mathematical community to accept such a "nice and rig- 
orous" definition, especially taking into account that its alternative in terms 
of infinitesimals was commonly viewed as "certainly non-rigorous" (although 
practically efficient). We can detect at least three more obvious reasons for 
this amazing phenomena in the history of calculus: 

1. In the period "from Leibniz to Weierstrass" , not only infinitesimals, 
but also real numbers, did not have a rigorous mathematical foundation. 
Thus, although the e, 5-definition makes perfect sense even in the framework 
of the rational numbers, this definition is completely fruitless without the 
completeness of the reals. We can proudly declare now that the obstacles 
related to the completeness of the real numbers belong to the past. Sadly, 
we can not be so proud about points 2 and 3 below. 

2. The e, 5-definition of limit is shockingly complicated due to the in- 
volvement of three non-commuting quantifiers V, 3, V. In the formulation for 
existence of a limit the quantifiers become four: V, 3, V, V. As a result, a real 
analysis course resembles a collection of exercises in mathematical logic rather 
than a rigorous version of calculus. The gap between the elementary calculus 
and real analysis widens and some students understandably wonder whether 
these two branches of mathematics have anything in common. Considerable 
efforts have been made to present the e, 5-definition in a more digestible and 
human-like form mostly by using geometric language (L. Gillman, R.H. Mc- 
Dowell j3] and S. Lang |6J). It is almost a public secret, however, that most 
mathematicians think and do research in terms of infinitesimals and use the 
e, 5-definition of limits only to present the final version of their work in a 
socially acceptable form. 

While acknowledging the importance of the above two factors for the 
unusually slow and late acceptance of the e, (^-definition in the history of 
calculus, we would like to focus our attention on another unpleasant feature 
of the e, 5-definition : 

3. The e, 5-definition of limit does not give any clue as to how to calculate 
the limit value L. At least this is true in the framework of the class of all 
functions (all polynomials or all rational functions, etc.), that is, functions 
with arbitrarily large derivatives. And this is exactly the situation in a typical 
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first calculus course. Thus, we have to guess a reasonable value for L and 
then prove or disprove the truthfulness of our guess with the help of the e, 5- 
definition. If the graph of the function is known (or, for the contemporary 
mathematician, if it is already on our computer screen), then the value of 
L can be reasonably guessed. Notice, however, that in this case we hardly 
need the concept of limit. In a numerical analysis course we probably would 
restrict our discussion to the class of functions with bounded derivatives (and 
given bound) and try to localize L within a given interval. But how to find L 
in a calculus course for freshmen who, presumably, do not know either what 
"limit" is, nor what "derivative" is, let alone the class of "functions with 
bounded derivatives" ? 

We sometimes are tempted to evaluate / for finitely many points different 
from r and try to guess the value of L assuming that there exists some 
pattern in the behavior of the function (J. Stewart [9], p. 50-61). Suppose, 
for example, that f(r + 1(T 10 ) = 1.99999999999. Maybe the correct limit of 
/ (as x approaches r) is L = 2 ? This guess is rooted in the following two 
myths: 

Myth 1: 1CT 10 is a "small number". More generally, "there are numbers 
in R which are small and others which are large". For example, 10 10 is "cer- 
tainly a large number". Indeed, we never use "millimeters" to measure the 
distance between two cities, nor do we count our annual income in "cents" . 
By changing the units of measurement (to suit our convenience) we always 
try to stay away from numbers such as 10~ 10 or 10 10 . The usage of different 
units for measurement explains the origin of this myth which, of course, has 
nothing to do with mathematics. 

Myth 2 : This myth originates in our experience as high-school students. 
It says that : "The integers are more likely to be the correct answers than 
are the fractions". Thus, we somehow prefer to believe that L = 2, not 
L = 1.999999999998, is the correct answer. 

Unfortunately, the values of a function / at finitely many points do not 
determine uniquely the limit of the function. We have to use even stronger 
language : The evaluation of a given function at finitely many points (differ- 
ent from the limit point r) is completely and totally irrelevant to both the 
value of the limit L and to the concept of limit in general. Here is the precise 
negative statement: 

Lemma 1 (Do not Guess !): Let f : X — > R, X C R, be a real function and 
r G R be a non-trivial adherent (accumulation) point of X . Let P(xi, yi),i = 
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1,2, ...,m, be finitely many (distinct) points in the plane M 2 such that Xj G 
X, Xi r and also such that Xi = Xj implies = yj . Let L G be an 
arbitrarily chosen real number (or even ±00 ). Then there exists a polynomial 
f : X — > R (or, a rational function g ) such that f(xi) — y^i — 1,2, ...,m, 
and lim^j. f(x) = L (or, \im x ^ r g(x) = ±00, respectively). 

Proof: : Choose a polynomial / of degree m and solve the linear system of 
equations f(xi) — y^i — 1, 2, m, f(r) = L, for the coefficients in /. The 
system has always at least one solution and we have lim^.-^ f(x) = L, since 
/ is continuous at r. ▲ 

In Section 2 we intend to show that if we use an arbitrarily chosen non-zero 
infinitesimal dx (instead of the increment 10~ 10 in our earlier example), then 
the value f{r + dx) uniquely determines the limit L. In Section 3 we present 
several examples from calculus to demonstrate how our method works in 
practice. For those readers who remain skeptical toward the practical merits 
of our approach, we remind them that all inventions of what we call today 
"calculus" have historically been discovered by means of infinitesimals. The 
reader who is interested in the history of calculus (C.H. Edwards, Jr. [2] 
and A. Robinson [S], Chapter X) will certainly observe that our method 
for calculating limits - if applied to calculating derivatives - is similar to the 
original Leibniz-Euler infinitesimal method. This explains the phrase: "Back 
to Classics" in the title of our article. 

This article is written for calculus teachers who are looking for alternatives 
to the conventional methods for teaching limits. We shall try to keep the 
exposition at the level of rigor in which complex numbers are defined: as 
"expressions of the form x+iy" , where x and y are real numbers and i = a/— 1. 
A more advanced (but still accessible) introduction to Infinitesimal Calculus 
is presented in Section of this paper, where the reader will find precise 
definitions and complete proofs. 

2 Preliminaries: Standard Part Mapping and 
Hyperreal Numbers 

We introduce the concept of infinitesimal and study the basic properties of 
hyperreal (nonstandard) numbers. We also study the operation known as 
the standard part mapping which is, in a sense, an algebraic counterpart of 
the concept of limit applied to numbers rather than to functions. 
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Definition 1 (Infinitesimals, Finite and Infinitely Large Numbers): 

(i) A number dx is called infinitesimal, in symbols, dx~0,if\dx\<l/n 
for all n G N. If x — y is an infinitesimal, we say that x and y are infinitely 
close, in symbols, x m y. 

(ii) A number x is called finite if \x\ < n for some n G N. 

(iii) A number x is called infinitely large if n < \x\ for all n G N. 

(iv) If x and t/ ^ are finite numbers, then the numbers of the form 
x/y are called hyperreal (or nonstandard) numbers. The set of all hyperreal 
(nonstandard) numbers will be denoted by *R. 

We denote by 2(*R), J(*R) and £(*R) the sets of thy infinitesimal, 
finite and infinitely large numbers in *R ; respectively. 

It is clear that all real numbers are finite and zero is the only infinitesimal 
in R. Also R has no infinitely large numbers. The following rules follow 
directly from the above definition: 

Theorem 1 (Properties): 

(i) finite ± finite = finite, finite x finite = finite. 

(ii) infinitesimal ± infinitesimal = infinitesimal. 

(iii) infinitesimal x infinitesimal = infinitesimal. 

(iv) infinitesimal x real = infinitesimal and, more generally, infinitesimal 
x finite = infinitesimal. 

(v) positive infinitely large + positive infinitely large = positive infinitely 
large. 

(vi) positive infinitely large x positive infinitely large = positive infinitely 
large. 

(vii) 1 /non-zero infinitesimal = infinitely large number. 

Warning: Numbers of the forms: ? n ^ es : ma ! , f^, gjjgffi ? argG , "positive 

° infinitesimal ' finite ' infinitely large ' r 

(negative) infinitely large - positive (negative) infinitely large" are always well 
defined provided that the denominators (if any) are non-zero. However, they 
can be of any type: infinitesimal, finite (real) or infinitely large. For example, 
let dx be a non-zero infinitesimal, in symbols, dx ^ 0, dx ~ 0. Then dx 2 /dx 
is infinitesimal, dx/dx is a real number (and that is 1), dx/dx 2 is infinitely 
large, both (2 + dx)/(2 + dx) 2 and (2 + dx) 2 /{2 + dx) are finite, dx~ l /dx~ 2 
is infinitesimal and dx~ 2 /dx~ l is infinitely large. Finally, 1/dx — 1/dx 2 is 
infinitely large, (1/dx + dx) — 1/dx = dx is infinitesimal and (1/dx + 5) — 
1/dx = 5 is finite (actually, real). 
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Remark 1 The level of rigor of Definition [TJ is similar to the level of rigor of 
the definition of complex numbers as "expressions of the form x + iy, where x 
and y are real numbers and i = y— 1 " . It is clear that our definition "hangs 
on" the existence of a non-zero infinitesimal (just as the existence of complex 
numbers "hangs on" the existence of % = y/—l). On the other hand, it is 
clear that if there exists one non-zero infinitesimal dx, then there are infinitely 
many infinitesimals: 2dx, 3dx, 4dx, dx 2 , dx 3 , etc. are also infinitesimals. 

Axiom 1: There exists a non-zero infinitesimal dx, in symbols, dx ^ 0, dx ~ 
0. 

Remark 2 The above definition together with Axiom 1 is a folk-like version 
of the statement: Let *R be a proper totally ordered field extension 

of R. We are simply trying to avoid fancy terminology. Recall that every to- 
tally ordered proper field extension of R is a non-Archimedian field, hence, it 
contains non-zero infinitesimals and infinitely large numbers. Strictly speak- 
ing, not any proper totally ordered field extension *R of R is adequate for 
the needs of Calculus. We also need that *R is a "non-standard extension" 
of R which means that *R satisfies two additional axioms. The first ax- 
iom (Transfer Principle) says, roughly speaking, that every function / in R d 
(d € N) has an extension */ in *R d such that the mapping * preserves the 
equivalence between the equations and inequalities in R and *R, respec- 
tively (where the right and left hand sides of these equations and inequalities 
are considered as real functions and their *-extensions, respectively). For 
example, we have —1 < sin a; < 1 •<=>■ x = x in R. Hence, by the Transfer 
Principle, it follows —1 < *sinx < 1 -<=>■ x = x in *R. In other words, the 
range of *sinx is the set { y G *R | — 1 < y < 1 }. We should mention that 
the concept of "nonstandard extension of a field" is in sharp contrast to the 
concept of an "algebraic extension of a field" (where the equivalence between 
some polynomial equations is, by design, violated in the extended field). An 
important consequence of the Transfer Principle is that *R is a real closed 
field, meaning that every polynomial equation of odd degree with coefficients 
in *R has a solution in *R. The second axiom (Saturation Principle) is a 
sort of completeness which implies, in particular, that every nested sequence 
of open intervals in *R has a non-empty intersection. For a more detailed 
exposition of nonstandard analysis by means of these two Principles we refer 
to the Keisler's Calculus textbook and its companion, written for calculus 
instructors (H. J. Keisler |l]-[5]). Both axioms (especially the Transfer Prin- 
ciple) are so natural that only an experienced mathematician will realize that 
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they are, actually, needed in Calculus. The situation is similar to the role of 
the Axiom of Choice in Real Analysis. It is well known that Real Analysis 
can not survive without the Axiom of Choice but it is completely possible to 
teach Real Analysis without even mentioning this axiom. 

If the reader still does not feel comfortable with the definition of hyperreal 
numbers, presented above, he/she should refer (now or later) to Section Oat 
the end of this paper. We resume our discussion on teaching calculus. It is 
clear that if r is a real number and dx is an infinitesimal, then r + dx is a 
finite number. Due to the completeness of R, the reverse is also true: 



Theorem 2 (Asymptotic Expansion): 
asymptotic expansion of the form x = 
and some infinitesimal dx ~ 0. 



Every finite number x in *R has an 
r + dx for some real number r 6 1 



Proof: See Theorem O in the last section of this paper. A 

Remark 3 (Completeness of R): H. J. Keisler [5], p. 17-18) proved that the 
statement of the above theorem is, actually, equivalent to the order com- 
pleteness of R (if R is treated merely as a "totally ordered field"). It is 
worth noticing that the completeness of the real numbers in the form pre- 
sented above appeared (treated as an obvious fact) in the early Leibniz- Euler 
Infinitesimal Calculus - 150 years before Cauchy, Bolzano, Weierstrass and 
Cantor formulated the completeness of the reals in the forms known from 
the contemporary real analysis textbooks. What does all this mean ? Well, 
perhaps the theory of the real numbers at the time of Leibniz and Euler was 
not so non-rigorous after all; it only takes so long until we finally figure out 
how the creators of Calculus preferred to express the completeness of R. 

In addition to the above properties of the finite numbers we have the 
following uniqueness result: 

Lemma 2 (Uniqueness): Let r 6 R and dx ^ 0. Then r + dx = implies 
both r = and dx = 0. 

Proof: r + dx = iff r = —dx. Hence, r = 0, since the zero is the 
only infinitesimal in R. ▲ 

The above property justifies the following definition. 
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Definition 2 (Standard Part Mapping): We define st : *R ->■ R U {±00} 

(a) st(r + cte) = r /or r£ R, drr 0; 

(b) st(l/dx) = ±00, dx ~ 0, /or drr > and drc < 0, respectively. 

In the case of finite numbers, the above several results can be summarized 
in the following corollary: 

Corollary 1 (Asymptotic Expansion of Finite Numbers): Every finite num- 
ber x G J-(*R) can be presented uniquely in the form x = st(x) + dx, where 
dx = x—st(x) is infinitesimal. We shall sometimes refer to the above formula 
as an asymptotic expansion of x. 

The next result follows immediately from the definition of "infinitesimal". 

Theorem 3 (Properties of st): Let x and y be finite numbers. Then we 
have: 

(i) x ~ y iff st(x) = st(y). In particular, st(dx) = for all in- 
finitesimals dx. 

(ii) st(r) = r for all r G R. 

(iii) Let x and y be not infinitely close. Then x < y iff st(x) < st(y). 
Consequently, for arbitrary x and y, "x < y implies st(x) < st(y)" and 
"x < y implies st(x) < st(y)". 

(iv) st(x ± y) = st(x) ± st(y). 

(v) st(xy) = st(x) st(y). 

(vi) st(x/y) = st(x)/ st(y) whenever st(y) 7^ 0. 

(vii) st(x n ) = (st( x)) n fo r all n e N. 

(viii) st(^/x) = \J st(x),n G N, where the condition x > (implying 
st(x) > ) is required in the case of even n. 

Remark 4 (Field, Ring, Ideal): In the usual algebraic terminology the re- 
sults of the above theorem can be summarized as follows: The set of hyperreal 
numbers *R is a totally ordered non-Archimedean real closed field, the set 
of finite numbers J r (*R) is a convex ring without zero-divisors (a totally 
ordered integral domain), the set of infinitesimals Z(*R) is a convex maxi- 
mal ideal in J r (*R) and the factor space J r (*R) F/X(*R) is isomorphic to 
R under st. 
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Remark 5 (Extended Real Line): If the standard part mapping st acts on 
infinitely large numbers, the result is either oo or — oo. In these cases we 
have to perform the usual legal and illegal operations in the extended real 
line R U {±00}: 

(a) Legal Operations: All Operations in K. are legal. In addition, the 
following are also legal (e is a positive real number). 

00 + 00 = 00, —00 — 00 = —00, ±e + 00 = 00, ±e — 00 = —00, 
00 x 00 = 00, (—00) x (—00) = 00, (—00) x 00 = —00, 
± e x 00 = ±00, e x (±00) = ±00, — e x (±00) = =Foo, 
1/ ± 00 = 0, ln(oo) = 00, e-°° = 0, e°° = 00. 

(b) Illegal Operations include: 

00-00, -, -, — , Ox (±00), 0°. 
' 0' 0' 00' V ; ' 

Here are several exercises with standard part: 

Example 1 (Exercises with st): In what follows r denotes a real number, 
and dx and dy are non-zero infinitesimals. 

1. st(dx) = 0. Similarly, st(3dx — Adx 2 ) = since 3dx — Adx 2 is an 
infinitesimal. 

r> lI -3+dx-dx 2 \ _ st(-3+dx-dx 2 ) _Q /O 
4. ^\2+2dxdy+dy A I st(2+2dx dy+dy 3 ) 

3 Incorrect- st ( yJ±<ibYjA = st( v / 3+^-y/3) = st(V3+di)- stjVs) = 

I dx J st(dx) st(dx) ' 

which is an illegal symbol. We disregard this calculation and try something 
else: Correct: 

y/3 + dx - y/3\ ( (V3 + dx - ^(^3 + dx + ^3) 

dx /V dx (V3~+dx + V3) 

( 3 + dx — 3 

= st 
= st 
= st 



Ax ( V3 + dx + V3) 
dx 

dx(V3 + dx + V3), 

1 \ strr 



.V^Tdx + VSJ st( v / 3 + dx + y/3) 
1 1 1 



^/st(3 + dx) + st( v / 3)) + 2V3 



9 



, which is the correct answer. 

4. Incorrect: st (^) = (^fgy = -3/0, which is an ille- 
gal symbol in R U {±00}. As before, we have to disregard this calcu- 
lation and try something else: Correct: st ( 2 ~dx+dx 2 ) = st ~2+dx ) = 
s ^ (ir) s ^ ( ~~2+dx ) = -*- 00 x ( — 3/2) = ±00, depending upon whether dx is 
positive or negative, respectively. 

Remark 6 (Guide): The following guide might help us to decide "what to 
do next" when we calculate the standard part st(x) of a hyperreal number 
x: 

(a) Let 1 be a finite number initially presented in the form x = r + dx 
(or it can be easily presented in this form). Then in order to calculate st(x) 
we have simply to apply the definition of st, i.e. to "drop the infinitesimal 
term dx" , st(x) = st(r + dx) = r. 

(b) Let 1 be a hyperreal number of unknown type (it might be a finite 
number but not presented in the form x = r + dx, or it might be an infinitely 
large number). To calculate st(x), we have to apply any of the properties 
of st, presented in Theorem [31 and/or any of the legal operations in the 
extended real line R U {±00} (Remark [5]). If the result of this calculation is 
a real number or a legal symbol in R U {±00} (part (a) in Remark [5]), then 
this is the correct answer for st(x). If at some stage of our calculations we 
obtain an illegal symbol in R U {±00} (part (b) in Remark [5]), we should 
stop, disregard the work done so far and start from the beginning a trying a 
different algebraic strategy. 

Remark 7 (Is the Algebra Familiar ?): The reader should not be surprised 
that the "algebra in the above examples sounds familiar" since both R and 
*R are real closed fields and, as we know, all real closed fields (Archimedean 
or not) obey the same laws of algebra. 



3 Limits Using Infinitesimals: Our Working 
Formulas 

In this section we derive several formulas for calculating the limit value L of 
a real function in terms of infinitesimals. Our working formulas - if consid- 
ered as definitions - are equivalent to the e, 5-definitions of the corresponding 
limits in real analysis but they do not involve a candidate L for the limit 
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value. Thus, we can start teaching a calculus course by "calculating" instead 
of "guessing and proving" . At the end of this section we summarize the main 
features of our method. 



Warning: The level of the following exposition is slightly higher than it is 
appropriate for teaching in class. To make it more accessible we recommend 
to following: 

(a) All details connected with the domain of the function should be 
skipped. They might be discussed later in the course when (and if) nec- 
essary. 

(b) The question of the existence of limits (which we discuss below) 
should be left aside or discussed when (and if) this question arises naturally 
in class. 

(c) The different types of limits should be presented in different sessions 
(not all at once as we have done below). 

(d) In the beginning the emphasis should be on those limits which have 
immediate geometric applications: the vertical and horizontal asymptotes. 
These limits are more important for the purpose of "sketching the curve" than 
the limits of the type lim x _;. r f(x) = L, where both r and L are real numbers 
(not ±00). The latter, although of fundamental importance for calculus, 
have a more subtle meaning, mostly to support the theory of "continuity" 
and "derivatives." 

Definition 3 (Limit): Let f : X — > R, X C R ; be a real function and 
(a, r) U (r, b) C X for some a,b,r G R,a < r < b. Suppose that 



for all non-zero infinitesimals dx and dy. Then st[/(r + dx)] G R U {±00} 
is called the limit of f as x approaches r, in symbols, 



where dx in the latter formula is an arbitrarily chosen non-zero infinitesimal. 
We shall refer to (TJ|) as our "working formula" for calculating limits through 
infinitesimals. 

Remark 8 (Existence of Limit): Notice that the condition (TjQ) guarantees 
the existence of the corresponding limit value st(/(r + dx)) = L and the 



(1) 



st[/(r + dx)] = st[/(r + dy)] in RU{±oo}, 



(2) 



lim/(x) = st[f(r + dx)], 
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independence of our working formulae on the choice of dx. We do not need 
to check the condition ([T]) before applying (J2J). Rather, we should start with 
02]) and when the calculations are done and the value L is obtained, we should 
check whether the value L depends on the choice of dx. If the answer is "no" , 
then L is the desired limit value. If the value of L depends on the choice of 
dx (say, L might depend on the sign of dx), then the limit lim x ^ r f(x) does 
not exist. As we already mentioned, it is preferable to skip the discussion of 
the existence of limit and focus on the working formula ([2]). 

Definition 4 (One Side Limits): 

(i) Let f : X — > R, X C R ; be a real function and (r, r + e) C X for 
some r G R and some e G R + and suppose that f has the property 

(3) st[/(r + dx)] = st[/(r + dy)} in R U {±oo}, 

for all positive infinitesimals dx and dy. Then st(/(r + dx)) G R U {±00} 
is called the limit of f as x approaches r from the right, in symbols, 

(4) lim f(x) = st[f(r + dx)], 

x— >r+ 

where in the above formula dx is an arbitrarily chosen positive infinitesimal. 

(ii) Let f : X — > R, X C R ; be a real function and (r — e, r) C X for 
some r G R and some e G R+. Suppose, in addition, that 

(5) st[/(r + dx)] = st[/(r + dy)] in RU{±oo}, 

for all negative infinitesimals dx,dy. Then st(/(r + dx)) G R U {±00} is 
called the limit of f as x approaches r from the left, in symbols, 

(6) lim f(x) = st(f(r + dx)), 

x—>r~ 

where dx in the last formula is an arbitrarily chosen negative infinitesimal. 

We refer to (fj]j and (0|) as our "working formulas" for the right and 
left-sided limits, respectively. 

Remark 9 (Left and Right Limits): The comparison between the above 
definitions implies that \im x ^ r f(x) exists iff each of lim a ._ ) . r + f(x) and 
linXj.^,,,- f(x) exists and lim a ._ >T + f{x) = lim^.^- f(x). In this case we have 

(7) lim f(x) = lim f(x) = lim f(x). 

x— >r x— >r+ x-» — 
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Definition 5 (Limit at Infinity): 

(i) Let f : X — > R, X C R, fre a rea/ function and (a, oo) C X /or some 
a G R. Suppose that 



(8) st(/(l/tte)) = st(/(l/dy)) zn R U {±00}, 

/or a// positive infinitesimals dx and dy. Then st(f(l/dx)) G R U {±00} is 
the limit of f as x goes to infinity, in symbols, 

(9) lim f(x) = st(f(l/dx)), 

x— >oo 

where dx is an arbitrarily chosen positive infinitesimal. 

(ii) Let f : X — > R, X C R, 6e a rea/ function and (—00,6) C X /or 
some b G R. Suppose that 

(10) st(/(l/tte)) = st(/(l/dy)) zn R U {±00}, 

/or a// negative infinitesimals dx and dy. Then st(f(l/dx)) G R U {±00} is 
the limit of f as x goes to minus infinity, in symbols, 

(11) lim f(x) = st(/(l/dx)), 

x— >— 00 

where dx is an arbitrarily chosen negative infinitesimal. 

We refer to (OJ) and ( TO)) as onr "working formulas" for the limits at 
infinity, respectively. 

Remark 10 (Existence of Limit at Infinity): As in the case of the usual 
limit (Remark |H]), the condition (JS]) or (JTP|) guarantees the existence of the 
limit value st(f(l/dx)) = L and the independence of the result on the choice 
of dx. As before we should start with the calculation of L by (JUJ) or ffTTj) . 
respectively, and when the calculations are over, we should check whether the 
value L depends on the choice of dx. If, not, then L is the correct answer. If 
L depends on the choice of dx, then the corresponding limit \im x ^± OQ f(x) 
does not exist. 

Remark 11 (Proper or Improper): If st(f(r + dx)) or st(/(l/dx)) is a real 
number, we say that the corresponding limits are proper. Otherwise (when 
a limit is 00 or —00), we say that the limit is improper. 
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Remark 12 (Unification): We shall often unite the right and left limits in 
the working formulae: 

(12) lim f(x) = st(/(r + dx)), lim f(x) = st(/(l/dz)), 

where in both formulae dx is an arbitrarily chosen infinitesimal, positive or 
negative depending on the sign in or in ±00, respectively. 

Theorem 4 (A. Robinson): The above definitions of different type of limits 
are equivalent to the corresponding e, 5 -definitions. 

Proof: : We refer the reader to Section [5] at the end of this paper. A 

4 Exercises on Limits 

The main advantage of the formulas st(/(r + dx)) and st(f(l/dx))) over 
the standard e, (^-definitions of lim^^i f(x) and lim^-too f(x), respectively, 
is that st(/(r + dx)) and st(f(l/dx)) prescribe an algorithm for calculating 
the limits in terms of / and r only (without involvement of a candidate L 
for the limit value) : 

(a) Evaluate / at the point r + dx (or at the point 1/dx, respectively), 
where dx is an infinitesimal, positive or negative, depending on the sign in 
r or in ±00, respectively. 

(b) Calculate the standard part of f(r + dx) (or the standard part of 
f(l/dx), respectively), following the rules in Section 1. 

Let us write once again our working formulae: 

(13) Iim/(z) = st[f{r + dx)\, 

x—>r 

where dx is an arbitrary non-zero infinitesimal, and 

(14) lim f(x) = st[/(r + dx)}, 

(15) lim f{x) = st[f{l/dx)], 

X— >±oo 

where in the last two formulae dx is an arbitrary infinitesimal, positive or 
negative, depending on the sign in r or in ±00, respectively. 

Here are several examples of applications of our working formulae. The 
reader will observe that our method requires less sophistication in factoring 
and less dependence on the Squeeze Theorem. 
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Example 2 \im x ^ r x n = st[/(r + cfo)] = st[(r + dx) n ] = (st(r + dx)) n = r n . 

Example 3 lim^ ±00 (l/a; n ) = st[f(l/dx)] = st{dx n ) = (st{dx)) n = 0. 

Example 4 \im x ^ ±(\x\/x) = st[f(dx)] = st(\dx\/dx) = ±1, where dx is 
an arbitrary infinitesimal, positive or negative, respectively. Notice that 
lim^olM/^) does not exist since lim^o+QxI/rr) ^ \rni x ^ Q - {\x\ / x) . 

Example 5 limbic*, (|rr| /a;) = st[f(l/dx)} = st(dx/\dx\) = ±1 (where, 
again, dx is an arbitrary infinitesimal, positive or negative, respectively). 

Example 6 

, im l£±i = st|/( -i +dx)]= st ^ 4 (- 1 + ^) + 1 



x^{-i)± x + 1 \ -1 + dx + l 

/-3 + Adx\ , , , . 
= st ( — J = st(l/dx) x st(-3 + 4cfe) 

= (±oo) x (-3) = =foo, 

where dx is an infinitesimal, positive or negative, respectively. 

Example 7 

,. x 3 + Ax 2 + x — 6 r „.„ , ., 
lim = st[/(l + dx)] 



x-*l X — 1 



st 



st 



4 + dx) 3 + 4(1 + dx) 2 + l + dx-Q 
1 + dx-l 
( 12dx + 7dx 2 + dx 3s 



\ dx 
= st(12 + 7dx + dx 2 ) = 12 

, where dx is a non-zero infinitesimal. These calculations will appeal to those 
students who prefer to expand the expressions (1 + dx) 3 and (1 + dx) 2 (and 
collect the like-terms), rather than to factor the cubic function x 3 +4x 2 — x— 6. 
Notice that the above calculations not only produce the correct limit value 12, 
but also present a rigorous proof that 12 is, in fact, the limit of the function, 
since the final result does not depend on the choice of the infinitesimal dx. 

Here is another example of an improper one-sided limit: 
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Example 8 



lim ( - X l I -) = st[f(2 + dx)] 

^2±\x 3 + x 2 -7x + 2J lJy n 

(2 + dxf - 9 



= st 
= st 
= st 



(2 + dxf + (2 + dxf - 7(2 + dx) + 2 
-1 + 12cfa; + 6cfo; 2 + ofe 3 

9<ix + 7dx 2 + dx 3 
-1 + 12dx + Qdx 2 + dx 3 



dx(9 + 7dx + dx 2 ) 



. 1\ f -1 + 12dx + Qdx 2 + dx 3 
— st — X st 



dx J \ 9 + 7dx + dx 2 
= (±oo) x (-1/9) = too 

, where dx is a non-zero infinitesimal, positive or negative, respectively. As 
in the previous example, these calculations will appeal to those students who 
prefer to expand the expressions (2 + dxf and (2 + dxf, rather than to factor 
the cubic function x 3 + x 2 — 7x + 2. 



Example 9 

x5±oo \ + 1 



lim ( ^F=^ ) = st[/(l/(fa)] 



= st 



1/dx 



{/l/dx 4 + l 

f\dx\ 1 
st j — - x 



dx yT+dx* 

\dx\\ ( 1 
1 X st ' 




1 



{/st(l + cfe 4 ) 
(±1) x 1 = ±1 

, where dx is an infinitesimal, positive or negative, respectively (notice that 
we use the fact that \/dx A = \dx\). 
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Example 10 lim !B _ > .± 00 [sin(a;)/a;] = st[f(l/dx)] = st(dx) x sm(l/dx)] = 
st(dx) x st[sin(l/<ix)] = x st[sin(l/dx) = 0, since st[sin(l/dx)] is a well 
denned (although explicitly unknown) real number in [-1, 1] (see Remark [2] 
about the range of sinx). 

Remark 13 (No Need of the Squeeze Theorem): We believe that the cal- 
culations in the last example, based on the product rule for standard part 
st (Theorem [2D, are shorter and simpler than the usual Squeeze Theorem 
arguments. Recall that the product formula for limits is non-applicable 
in the case of the function sinx/x since lim^oo sin x does not exist in 
K. (even as an improper limit). In contrast, we have st[dx sin(l/ofcc)] = 
st(dx) st[sin(l/dx)] since each of st(dx) and st[sin(l/cfe)] exists in R, by 
Theorem [2] and Lemma [2J In general, the infinitesimal method for calculat- 
ing limits is less dependent on the Squeeze Theorem, and more rarely requires 
the use of inequalities, compared with the usual Weierstrass e, 5-method. We 
consider that feature of the Infinitesimal Calculus as an advantage over the 
conventional standard methods, taking into account that the students in 
calculus are rarely in love with inequalities. 

Although our text is devoted to limits only, we shall shortly mention the 
definitions of derivative and differential in terms of infinitesimals: 

Remark 14 (Derivative): We define the derivative by 

= st (f±±M^r 

\ (JjJb 

where dx is an arbitrary non-zero infinitesimal. For example, 

x + dx) 3 — x 3 \ fx 3 + 3x 2 dx + 3xdx 2 + dx 3 — x 3 



(* 3 ) = st " iz = st 



st 



dx J V dx 

'3x 2 dx + 3xdx 2 + dx 3 



dx 



- ) = st (3x 2 + 3xdx + dx 2 ) = 3x 2 . 



Similarly, we define differential by dy = f'(x)dx, where dx is an infinitesimal. 
Thus, the Leibniz notation dy/dx = f'(x) holds "by the definition" of dy for 
all non-zero infinitesimals dx. 
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Summary: We summarize the properties of our working formulae (jl3p - (fT5|) : 

1. Our working formulae, if considered as definitions, are equivalent of 
the usual e, ^-definitions of the corresponding limits at any level of generality 
and rigor (see the next section). 

2. Our working formulae are free of a candidate, L, for the limit value. 
Hence we do not need to guess and check (because there is nothing to guess). 
Also, our working formulae "really work" in the sense that they produces the 
correct value of L, as demonstrated by the above examples. 

3. Our method requires less sophistication in factoring (Examples flTJ)- 
flB]) and it is less dependent on the Squeeze Theorem compared with the 
conventional methods (Example (ITU]) and Remark (1131) ). 

4. Under the assumption that "the limit exists", our working formu- 
lae (fT3"]) - ([T5]) are free of quantifiers, as opposed to the three non-commuting 
quantifiers "V, 3,V" in the e, 5- definition of limit. On the other hand, each 
of the criterions for existence of limit ((ITJ), fl3]), ()5]), OH]), ( ITO]) ) involves two 
commuting quantifiers "V, V" only, as opposed to the four non-commuting 
quantifiers "V, 3, V, V" in its standard counterpart. Thus, our method is eas- 
ier to apply to rigorous proofs when (and if) the teacher decides to do rigorous 
proofs. In fact, in our method, the rigorous proof that "L is, actually, the 
correct limit value" coincides with the "calculation of L." 

As a result, the Calculus becomes to "calculate" again, as it was originally 
designed to do. 

Remark 15 (Infinitesimals in Mathematica): Assume that you are already 
"addicted" to calculators and computers and that you are not planning to 
"quit" any time soon. Assume that your students have already purchased ex- 
pensive calculators and are eager to calculate "anything which comes along." 
Assume that your university has already spent a lot of money on buying 
computers and the spending has to be somehow justified. And to complete 
the scenario assume, finally, that you have just won a generous grant from 
NSF for "using technology in teaching calculus." Under these circumstances 
your dilemma will be "how to reconcile the computers with infinitesimals"? 
The good news is that computers are able to handle infinitesimals, and they 
actually work with infinitesimals anyway. Take, for example, "Mathemat- 
ica". Have you ever thought about how Mathematica calculates limits ? It 
might occur to you that the computer evaluates the function / at finitely 
many points and announces one these values for the "correct answer" ? 
Or, perhaps, the computer has simply memorized the limits of all possible 
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functions "you will ever ask it for"? The answer is "neither of the above." 
Rather, Mathematica calculates the formal (Taylor) asymptotic expansion of 
f(r + dx) by the command "Series", treating dx as a "formal variable" and 
truncates all terms in the series but the first by the command "/. dx — > 0". 
This procedure has very little to do with the e, 5-definition of limit and it 
is almost identical to the operation "taking the standard part" discussed in 
this article. The framework of these calculations is the field M.((dx)) of for- 
mal Laurent series with real coefficients and formal variable, denoted by dx. 
Notice that the field is non-Archimedian and the formal variable, 

dx, if considered as an element of R((cfc)), is a positive infinitesimal. So 
what ? Well, it means that Mathematica (believe it or not) calculates limits 
through infinitesimals in the framework of the field M.((dx)). The formula 
st[f(x + dx)] = Series[f(x + dx), dx, 0, 1}/ .dx — >■ can be used for calculat- 
ing the "standard part" of f(x+ dx) in Mathematica if you decide to do so. 
The author of this article, however, is unable to see any pedagogical merits 
of this activity unless, perhaps, for the purpose of a better understanding of 
how Mathematica works. 



5 An Introduction to Infinitesimal Calculus 

Here we present a short introduction to the modern Infinitesimal Calculus 
known as well as A. Robinson's Nonstandard Analysis. We would like to 
assure the reader that the usual background in real analysis is more than 
enough to follow this text. For more detailed exposition we shall refer to 
(H. J. Keisler |l]-[5]) and (T. Lindstr0m [7]), where the reader will find more 
references to the subject. For a really short (although somewhat dense) expo- 
sition of both axiomatic and sequential approaches to nonstandard analysis 
we refer to T. Todorov [10], p. 685-688. We shall restrict our exposition to 
the nonstandard treatment of proper limits of the form lim;^,. f(x) only, and 
leave the improper limits, as well as the limits at infinity, to the reader. All 
results in Section 2, follow as particular cases. 

Although the nonstandard analysis arose historically in a close connection 
with model theory and mathematical logic (A. Robinson [8]), it is completely 
possible to construct it in the framework of the standard analysis, i.e. as- 
suming only the properties of the real numbers (along with the Axiom of 
Choice). The method (due to W. A. J. Luxemburg), is known as the ultra- 
power construction or constructive nonstandard analysis: 
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1. Let N be the set of the natural numbers and V( N) be the power set 
of N. Let \i : P(N) — > {0, 1} be a two-valued finitely additive measure such 
that n(A) = for all finite A C N and = 1. We shall keep /i fixed in 
what follows. 

Remark 16 (Existence of /i): To show that there exists a measure with 
these properties, it suffices to take a free ultrafilter Won N and define \i by 
fJ>(A) = 1 for A E hi and fi(A) = for A hi. Recall that a non-empty set 
hi of subsets of N is called a free ultrafilter on N if it satisfies the following 
four properties: (a) hi is closed under intersection; (b) li A,B C N, then 
hi 3 A C B implies B Ehi; (c) For any A C N exactly one of the following 
is true: A e hi or N \ A e hi; (d) f] AeU A = 0. Recall that the existence of 
free ultrafilters on N, follows from the Axiom of Choice (H. J. Keisler [5], p. 
49). We should mention that the familiarity with the theory of ultrafilters is 
not necessary for the understanding of what follows. 

The next properties of fi follow immediately from the definition: 

Lemma 3 (Properties of fi): Let A, Be N. Then: 

(a) fi(A U B) = 1 -<=>■ [/J>(A) = 1 or fi(B) = 1 ]. In particular, for any 
A C N exactly one of n(A) = 1 and //( N \ A) = 1 is true. 

(b) n(A) = 1 for all co-finite sets A of N. In particular, /i(N) = 1. 

(c) fi(A UB) = lfi(A) = and fx(B) = 0]. 

(d) n(A) =n(B) = 1 v{AnB) = 1. 

(e) A C 5 C N and fi(A) = 1 implies fi(B) = 1. 

2. Let M N be the set of all sequences of real numbers considered as a ring 
under the usual pointwise operations. Define an equivalence relation ~in 1 N 
by: (a n ) ~ (b n ) if a n = b n a.e. (where "a.e." stands for "almost everywhere"), 
i.e. if fi({n \ a n = b n }) = 1. Then the factor space *M = M N / ~ 
defines a set of nonstandard real numbers (or hyperreals). We shall denote 
by (a n ) the equivalence class determined by the sequence (a n ). We also 
define the embedding R C *M by r — >• (r, r, r, . . . } . In what follows we 
shall identify notationally a given real numbers r with its image in *R. 
The addition and multiplication in *M is inherited from M N . The order 
relation in *M is defined by: (a n ) < (b n ) if a n < b n holds a.e., i.e. if 
^({n | a n < b n }) = 1. We define also |x| = max {x, — x}. Notice that we 
define one specific nonstandard extension *M of the reals R which depends, 
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in general, on the choice of the measure fi. We should mention that the 
different fields of the form *R (corresponding to different measures /x) are 
not necessarily isomorphic to each other. We also have card(*R) = card(R). 

Theorem 5 *R is a totally ordered non- Archimedean field containing R as 
a totally ordered subfield. 

Proof: *R is a ring since R N is a ring. To show that *R has no zero 
divisors, assume that (a n ){b n ) = in *R, i.e. /i({n | a n b n = 0}) = 1. 
Denote A = {n \ a n = } and B = {n \ b n = 0} and observe that 
{n | a n b n = 0} = AU B, since R has no zero divisors. Hence either 
fJ>{A) = 1 or /i(-B) = 1, by Lemma [3J i.e. either (a n ) = 0, or (b n ) = 0, 
as required. To show that the non-zero elements in *R are multiplicative 
invertible, assume that (a n ) ^ in *R, i.e. /x({n | a n 7^ 0}) = 1. Denote 
{n I a n 7^ 0} = C and define (6 n ) G R N by b n = l/a n if n G C and anyhow 
(say, 6 n = 1) if n G N \ C. We have C C { n |a n 6 n = l} which implies 
/i({n I a n 6 n = 1 }) = 1, by Lemma [3J Thus (a n )(b n ) = 1, as required. 
To show the trichotomy of the order relation, assume that (a n ) 7^ (b n ) and 
denote A = {n \ a n < b n }. We have N\A={n \ a n > b n } and, thus, 
exactly one of fi(A) = 1 or /x(N \ A) = 1 is true, by Lemma [3j That is 
(a„) < (b n ) or (a n ) > (b n ) which is equivalent to (a n ) < (b n ) or (a n ) > (b n ), 
as required, since (a n ) 7^ (6 n ), by assumption. The rest of the properties of 
the totally ordered field can be proved similarly. The embedding R C *R is 
obviously field and order preserving. To show that *R is non- Archimedean, 
observe that m < (n) in *R for any m in N, (where m is considered as an 
element of *R) since the set { n \ m < n } is co-finite, and hence, of measure 
1. ▲ 

Example 11 (1/n), (1/n 2 ), (1/ Inn), (e~ n ) are positive infinitesimals (differ- 
ent from each other) and (n), (n 2 ), (Inn), (e n ) are positive infinitely large 
numbers (also different from each other). The number (3 + 1/n) = 3+ 
is finite (but not real). Let us take the first example: for any m G N, the 
set {n I < 1/n < l/m} is co-finite, hence, of measure 1. Therefore, 
< (1/n) < 1/m in *R, i.e. (1/n) is a positive infinitesimal. The rest of 
the examples are treated similarly. 

3. It is clear that R C J 7 (*R),X(*R) C J r (*R), RnZ(*R) = {0}, J"(*R)n 
£(*R) = and J 7 (*R) U £(*R) = *R. From the above definition it fol- 
lows easily that J r (*R) is a totally ordered integral domain and X(*R) is 
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a convex maximal ideal in J r (*IR). Hence, J r (*R)/X(*IR) is a totally or- 
dered field which is isomorphic to R as totally ordered fields. The canon- 
ical homomorphism st : J-"(*R) — > R is called the standard part map- 
ping. Notice that st(a n ) exists for any bounded sequence (a n ) in f N and 
st(a„) = lim^oo a n for any convergent (a n ). Conversely, if (a n ) is a finite 
number, then st(a n ) = lim^oo for some subsequence (ak„) of (a n ) such 
that k n | n G N}) = 1. The following result follows immediately: 

Theorem 6 (i) Let x G *R. Then x G J r (*M) iff x = r + dx for some 
x G R and some dx G Z(*IR). 

(ii) If x E J r (*R) ; then the presentation x = r + dx is unique and r = 
st(x). In particular, st(r) = r for any re R. 

(iii) The standard part mapping is order preserving in the sense that 
x < y in .F(*R) implies st(x) < st(y) in R. 

4. Let X C R. Then the set *X = { (x n ) G *R | x n G X a. e. } is called 
the nonstandard extension of X. For any X C R we have X C *X and 
X = *X iff X is a finite set. The above definition holds also in the case 
when X C R d (d G N). If X C R and F C R, then *(X x F) = *X x *F. 
In particular, we have *(R d ) = (*R) d , so we can write simply *R d . 

Example 12 Let Q, Z, N, E, O, P, etc. be the sets of the rational, integer, 
natural, even, odd, prime, etc. numbers, respectively. Then the elements 
of *Q, *Z, *N, *E, *0, *P, etc. will be called nonstandard rational numbers, 
nonstandard integer numbers (hyperintegers), nonstandard natural numbers 
(hypernatural numbers), etc., respectively. The set of the infinitely large 
natural numbers, i.e. the infinitely large numbers in *N, will be denoted by 
Noo. We have *N = N U Noo. If a, b G R, then 

*{xe R \ a <x <b} = {x e *R |a<a;<6}, 

which will be denoted for short by *[a, b]. It follows that r + dx G *[a,b] for 
all r G [a, b) and all non-negative infinitesimal dx. Similarly, we have 

*{x e R \ a < x < b} = {x e *R | a < x < b}, 

which will be denoted for short by *(a, b). We have r + dx G *(a, b) for all 
r G (a, 6) and all infinitesimal dx . 
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Theorem 7 (Adherent Point): Let r G R and X C R. Tnen r is a non- 
trivial adherent point of X if f there exists dx G *R sncn tnai dx 7^ 0, dx ~ 
and r + dx G *X (or, equivalently, iff there exists x G *X snc/i iaai r/i 
and r m x). 

Proof: : (=>) For any n G N the set X n = { x G X | < \x - r\ < ^ } 
is non-empty, by assumption. Hence, by the Axiom of Choice, there exists 
(x ra ) in R N such that x n G X n for all n G N. Now, dx = (x n ) — r is the 
infinitesimal we are looking for. Indeed, we have < \dx\ < (-), hence 
dx 7^ and dx m 0, since (-) ~ 0. Also r + dx = (x n ) G *X, since x n G X n 
for all n G N. (<=) We have dx 7^ 0, dx « 0, r + dx G *X for some dx G *R, 
by assumption. Suppose m G N and observe that < \dx\ < 1/m. We have 
dx = (e n ) for some (e n ) in R N . The set { n | < |e n | < 1/m, r + e n G X } is 
of measure 1, hence, it is non-empty. The latter means that r is a non-trivial 
adherent point of X. k 

5. Let / : X — > R be a real function, where X C R. Then the function 
*f : *X ->■ *R, defined by */((x n )) = (/(x n )) for all (x n ) G *X, is called 
the nonstandard extension of / since */(r) = /(r) for all r G X. The above 
definition holds also in the case X C R d (d G N). 

Theorem 8 (A. Robinson): Let r be a non-trivial adherent point of 
f : X — > R, X C R ; be a real function and L G R. Taen Hm a; _ >r /(x) = 
2// *f{ r + dx) ~ L for all dx G *R such that dx 7^ 0, dx ~ and 
r + dx G *X. // t/ie /zmzt exists in R ; iaen lim^j. /(x) = st(*/(r + dx)). 

Proof: : (=>) Let e G R+. By assumption, there exists 5 G i?+ such that for 
all x G X, < |x — r\ < 5 implies \ f{x) — L\ < e. Let dx 7^ 0, dx ~ and 
r + dx G *X for some dx G *R. Notice that dx exists, by Theorem [3, since 
r is a non-trivial adherent point of X, by assumption. We have dx = (e„) 
for some sequence (e n ) in R N . Next, we define the sets: 

A s = {n I < \e n \ < 5 and r + e„ G X } and B € = {n \ \f(r + e n ) — L\ < e }. 

We have fi(B e ) = 1, by Lemma El since As C £> e and = 1, by assump- 

tion. Recapitulating, we have \*f(r + dx) — L\ < e for all e G R+, which 
means that *f(r + dx) ~ L, as required. 

(<=) Assume (on the contrary) that lim x ^ r f(x) = L is false. Thus, there 
exists e G R+ such that the sets 

X n = { x G X I < \x — r\ < - and \f(x) - L\ > e }, 
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are non-empty for each n G N. By the Axiom of Choice, there exists a 
sequence (x n ) in R N such that x n G X n for all n G N. Define dx G *R 
by = (x n ) — r. We have < \dx\ < (-) hence dx 7^ and da; ~ 0. 
Also r + dx = (x n ) G *X and |*/(r + dx) — L\ > e. The latter means 
that *f(r + dx) — L is a non-infinitesimal, a contradiction. The formula 
\im. x ^ r .f(x) = st(*/(r + da;)) follows directly from *f(r + da;) ~ L after 
applying the standard part mapping to both sides and taking into account 
that st(L) = L since L is a real number. A 

In order to eliminate completely the "candidate" for limit L from our 
theory, we have to present a nonstandard characterization of the existence of 
a proper limit in terms of infinitesimals similar to the Cauchy convergence 
criterion: 

Theorem 9 (Existence): Let f : X — > R, X C R ; be a real function and 
let r G R be a non-trivial adherent point of X . Then the following are 
equivalent: 

(i) The limit lima,-;.,, f(x) exists in R. 

(ii) / is fundamental (or Cauchy) toward r in the sense that 

(Wee R + )(3<5g R+)(Vx,y G X)[0 < \x-r\, \y-r\ < 6 \f(x) - f(y)\ < e). 

(iii) *f(x) ~ *g(y) for all x,y G *X such that x 7^ r,y 7^ r,x ~ r and 
y « r. 

(iv) */(r + dx) ~ *p(r + dy) /or all dx,dy 7^ 0,dx,dy ~ suc/i t/iat 
r + dx, r + dy e *X. 

(v) st(*/(r + dx)) = st(*/(r + dy)) G R (Tmt never become ±00) for all 
dx, dy 7^ 0, dx, dy ~ such that r + dx, r + dy G *A. 

(vi) (3/i G Z(*M + ))(Vdx,dy G *R) 

[r + dx,r + dy e *X and < |dx|, \dy\ < h] [*/(r + dx) ps */(r + dy)] , 

where I(*W + ) denotes the set of the positive infinitesimals in *R. 

Proof: : (i) <^ (ii) is the Cauchy Criterion for existence of limits (Alan F. 
Beardon p], Theorem 4.4.1, p. 57). 

(i) =Kiii): We have *f(x) « L and */(?/) ~ L for the same L G R, by 
Theorem [HI hence, */(x) ~ *f(y), as required. 

(iii) -v4>(iv) follows immediately by letting x = r + dx and y — r + dy. 

(iv) => (vi) in a trivial way. 
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(vi)=^(i) : We have h = (h n ) for some sequence (h n ) in R N . Without loss 
of generality we can assume that h n > for all n G N. Now, suppose (for 
contradiction) that (i) fails, i.e. there exists e G R + such that A n ^ for all 
n G N, where 

4 = {(i^)GlxI < \x-r\, \y-r\ < h n and \ f(x) - f(y)\ > e}. 

Hence (by Axiom of Choice), there exists a sequence (x n ,y n ) in X N x X N 
such that (x n ,y n ) G A n for all n G N. We define the nonstandard numbers 
(x n ), (y n ) G *X and observe that 

0<|(x n )-r|, \(y n )-r\<(h n ) and |*/<z n » - *f«y„»i > e, 

in *R, by the choice of (x n ) and (y n ). Thus, *f((x n )) — *f((y n )) is a non- 
infinitesimal, which contradicts (vi) for ft, = (h n ), dx = (x n ) — r and dy = 
(Vn) - r. 

(i)^>(v) : We have st(*/(r + <ix)) = L and st(*/(r + dy)) = L for the 
same L G M, by Theorem [HI hence, st (*/(#)) = st (*/(?/)) e as required. 

(iv) =>■ (v) : st(*/(r + dx)) = st(*/(r + dy)) G R implies, in particular, 
that */(r + da;) and *f(r + dy) are finite numbers, thus, */(x) ~ */(y) follows. 
A 

Simplified Notation: For the purpose of teaching and explicit calculations 
we recommend the following simplified notations: 

(a) We shall skip the asterisks in front of */, writing simply /. This is 
perfectly justifiable since */ is an extension of /. 

(b) If X C R, then we shall sometimes write simply X meaning *X. For 
example, we shall write [a, b] meaning, actually, 

*[a,b] = {x e *R \a<x<b}. 

(c) Finally, we shall write (—00, 00) for both R and *M leaving the 
reader to figure out from the context which one we mean. 

(d) We prefer to use the terminology "hyperreal numbers" rather than 
"nonstandard numbers" (to avoid the shocking effect of the word "nonstan- 
dard"). 

(e) We preserve our rights to come back to the more precise *-notation 
when (and if) needed. 
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